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I. INTRODUCTION 



The concepts of phase transition and critical points are denned, strictly speaking for macroscopic systems. However, 
it has been recently suggested that precursors of phase transitions can be observed in finite-size mesoscopic systems . 
In Nuclear Physics the different nuclear shapes and the phase transitions between them are conveniently studied within 
the Interacting Boson Model (IBM) Q ■ This was early recognized after the introduction of the model 0, 0, IE IE 
but has been studied thoroughly in the last few years HIEllEIHIHEElllEIIEIHIH after the introduction 
of the concept of critical point symmetries [T^. l2fl l2l| . Since the IBM is formulated from the beginning in terms of 
creation and annihilation boson operators, its geometric interpretation in terms of shape variables is usually done by 
introducing a boson condensate with two shape parameters, (3 and 7 (order parameters) [EE2- The parameter (3 is 
related to the axial deformation of the system, while 7 measures the deviation from axial symmetry. The equilibrium 
shape of the system is obtained by minimizing the expectation value of the Hamiltonian in the intrinsic state. Shape 
phase transitions are studied theoretically using one (or a few) control parameter(s) in the Hamiltonian. These control 
parameters drive the system in different phases characterized by order parameters and allows one to study in a simple 
way phase transitions and critical points in Nuclear Physics. 

The phase diagram of the IBM has been studied with several approaches HIEIIEIllHQIIEIHm 

and it is 

well known that the dynamical symmetry associated to U(5) corresponds to a spherical shape ((3 — 0), the dynamical 
symmetry SU(3) is associated to an axially deformed shape (7 = 0, 7r/3, f3 7^ 0) and the dynamical symmetry 0(6) 
is related to a 7-unstable deformed shape ((3 7^ and 7-independent). These symmetry limits are usually represented 
as the vertices of a triangle (Casten triangle) • Phase transitions between these shapes have been widely studied 
and it is known that the phase transition from U{h) to 0(6) is second order while any other transition within the 
Casten triangle from a spherical to a deformed shape is first order [E0|- These studies have been performed, as 
mentioned above, by using the intrinsic state formalism. However, it is known that this approximate method is only 
correct at leading order in a 1/N expansion where N is the number of bosons. In this paper, we present a method 
to go beyond this order and compute finite-size corrections to several spectroscopic observables. We stress that 1/N 
corrections obtained with the intrinsic state formalism (or Hartree-Bose method) are in general incorrect and give no 
information on the proper finite-size corrections. 

The paper is organized as follows: first the model Hamiltonian is introduced in Sect. II. In Sect. Ill, the Holstein- 
Primakoff mapping |24| is performed leading to a boson Hamiltonian in which we retain terms in orders N, N 1 / 2 and 
TV . Then a Bogoliubov transformation is performed to diagonalize the Hamiltonian and to study both the symmetric 
(spherical) and the broken (deformed) phases. All this is done in general for two- level boson models in which the 
lowest level is a scalar (s) boson while the upper level is an arbitrary L boson. The IBM corresponds to the particular 
case L = 2 (d^ bosons). In addition to the IBM, we present results for the case L = as an illustration of the general 
method. In Sect. IV we compare the analytical results with exact numerical diagonalizations for different paths along 
the Casten triangle. Finally, Sect. V is for the summary and conclusions. 
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II. THE MODEL 

As it has been noted before [Llj , the experimental exploration of the shape transition and critical points in nuclei is 
difficult due to the lack of a continuous control parameter. However, in theoretical studies this limitation is overcome 
by using a Hamiltonian written in terms of one or more control parameters that can vary continuously. In this work, 
we consider a two-level boson model in which the lowest level is characterized by a zero angular momentum (s-boson) 
while the upper level has an arbitrary angular momentum L. The Hamiltonian proposed is a generalization of the 
IBM consistent-Q formalism (CQF) (23, which depends on two control parameters x and \ 

H = xn L -^QX-QX, (1) 

where = L^L^ is the operator for the number of bosons in the upper level, N is the total number of bosons, 

the symbol • stands for the scalar product defined as a ■ b — ^/t—_i J {— l)** /^-/*) an( l Q x ' IS a multipole operator 
written as, 

Q* = ^L+LU)f)+xpxL}£\ (2) 

where = (— 1)^L_^. For L = 2, (d-bosons) the Hamiltonian Q is the well-known CQF Hamiltonian for IBM. 
Though it is not the most general IBM Hamiltonian, it captures the most important low energy properties of a wide 
range of nuclei j26l I27I |28| . In particular, it is general enough to describe different nuclear phases and quantum phase 
transitions, and it has been used for that purpose at the mean field level |9llTollT2|. 

The Hamiltonian comprises different models depending on the value of L. For instance, for L = 1 the Hamilto- 
nian is appropriate for studying the phase diagram of the vibron model |29| of interest in Molecular Physics. 

III. MEAN FIELD AND BEYOND 

The usual way of getting the phase diagram of the model is to introduce shape variables. This can be done by 
considering the intrinsic state formalism, also called Hartree-Bose approximation, |3l l5l 1221 1301 ] . In this approach, the 
ground state is a variational state built out of a condensate of "dressed" bosons, that are independent bosons moving 
in the average nuclear field. For L = 2, these bosons are defined as 

' s f +/3cos7<4 + -L#siri7(4 + d - 2 )V ( 3 ) 



and the N boson condensate is, 



\c) = -L=(rtno). (4) 



The variational variables f3 and 7 are the order parameters of the system and their equilibrium values are fixed 
by minimizing the expectation value of the energy. The expression of this energy can be found in many references 
[3I I22I I30L l3l | and can be written schematically as follows, 

E(N, /?, 7, x, X ) = NF^(N, 0, 7, x, + l)F^(N, ft 7, x, x ), (5) 

where F^(N, /3, 7, x, \) ' IS the matrix element of the one-body operators divided by N and F^(N, 0, 7, x, \) is 
the matrix element of the two-body operators divided by AT — 1. Note that there is no iV 2 dependence in the two-body 
operator due to the definition of the Hamiltonian. Actually, the only relevant contribution is the leading one (order 
N) since the next one (N° for instance) are incomplete as explained below. 

For the standard IBM Hamiltonian (L = 2), with an attractive quadrupole interaction, the nucleus always becomes 
axially deformed, either prolate (7 = 0) for \ < or oblate (7 = tt/3) for x > 0. As a consequence, the parameter 
7 can be incorporated in the value of j3. f3 > corresponds to 7 = while negative (3 implies 7 = tt/3. In the case 
X = the nucleus becomes 7 unstable, i. e., the energy is independent of 7. 

In this framework, one-phonon excitations above the ground state are constructed by directly replacing in the 
ground state H) a condesate boson by an excited boson (TDA method) or by including ground state fluctuations 
(RPA method) |3. Isil l82l | . For L — 2, there are five excited phonons that are characterized by their angular momentum 
projection K and can be labelled as: /3-excitation with K = 0, 7-excitations with K — ±2 and finally two K = ±1 
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excitations. It should be noted that not all the excited phonons are always physical, some of them become spurious, 
Goldstone bosons associated with broken symmetries. This is the case for axially deformed nuclei, the K = ±1 
excitations are Goldstone, spurious bosons because the state constructed with this excitation corresponds to a 0(3) 
rotation of the whole system. In the case of 7 unstable nuclei, the K = ±2 excitations also become Goldstone bosons 
and are related with 0(5) rotations of the ground state. In the case of L — only a K — excitation exists and it is 
directly related, as we will see, with the band of the IBM 

The mean field description of the ground state energy just mentioned is only valid at order N. The first quantum 
corrections can be obtained within the RPA formalism. Alternatively, the Holstein-Primakoff expansion [24| offers a 
simple and natural expansion in powers of l/N. The advantages of this transformation are: it is Hermitian, preserves 
the boson commutation relation, provides a correct expansion in powers of N and its leading order coincides with the 
mean field contribution. 

The Holstein-Primakoff expansion consist in eliminating the s-boson transforming the bilinear boson operators in 
the following way, 



b\b v , 

^/2 6 t (1 _ n6/iV )V2 = (s t^ ) t ) 

N - n b . 



(6) 
(7) 
(8) 



where the 6-bosons satisfy [6^,6^,] 



The mapping fulfils the commutation relations at each order in N in the 



Taylor expansion of the square root. 

We next introduce the c-bosons through a shift transformation 



bl 



NX* 



(9) 



where the A^'s are complex numbers which form a (2L + l)-dimensional vector. This shift allows for a macroscopic 
occupation number ni,. Thus, it allows to consider at the same time the spherical, setting A M = for all /x, and the 
deformed phase, A^ 7^ 0. In this latter situation, we shall only consider the case Ao ^ without loss of generality. 
The Hamiltonian then reads, 



H = N 



^ {5x - 4 - 4(x - 1)X 2 + (x- 1)x45a [4(1 - X 2 Q f' 2 + X <*$>*] } + 



N^Xo (c + c ) 1 5a; - 4 - 8X 2 Q (x - 1) + 2(x - l)x^M [ (1 _, A ^~|i/2 + X a ifi X ° 
N°\ [3x - 2 - 6X 2 (x - 1)] n c + (x - 1) [(2L + 1) - (2L + 3)A§ + (l - Aq) (Pj + P c ) - 4.X 2 [c\ 2 + 2c c + c, 



2 X (x-l)\x (l-Xl) 1 / 2 



+ L 

£< 

fj,— — L 



(L) 
*0,P 



A a 0.0 




3 c 



2cjc 



2n, 



A a 0,0 



4(l-Ag)3/a 



LL— — L 

= H 1 +H 1/2 + H + O(l/^), 




where aj^l — (L, fx; Lv\L. fi + v) and Pj 
1 of order N ( 

E(N, (3, x,y)=N 



(Pc) 1 



The term of order N (Hi) is exactly the mean- field energy. Setting A = /3/y/l + (3 2 one gets, 

1 



(1 + /? 2 ) 5 



5a; 



-xp 2 + f3y(x- l)(A + (3y) 



(12) 



where y = X a oo- I n the case of L = 2, equation (|12|l reduces to the IBM ground state energy. Note that Hi only 
depends on L through the Clebsch-Gordan coefficient (L, //; Lv\L, /z + v) although this dependence can be absorbed 
in the parameter y. 
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spinodal critical antispinodal 



FIG. 1: Qualitative phase diagram for the Hamiltonian Q and L = 2. The insets show typical energy surfaces versus the 
deformation parameter f3 in each of the phases and at the phase borders. The control variables of the diagram are defined as 
p=l-a:and0 = ^= (^r+x). 



Hi provides the mean field energy and therefore the equilibrium values of the order parameter. In Fig.^is depicted 
the phase diagram corresponding to Hi. For given parameters x and x(y), the first step consists in minimizing Hi 
with respect to /3(A) getting the equilibrium value /?o(Ao). The study of these minima have been shown in several 
publications 0,0, but for completeness we summarize here the main features: 

• [3 = is always a stationary point. For x < 4/5, (3 — is a maximum while for x > 4/5 becomes a minimum. 
In the case of x = 4/5, (3 = is an inflection point, x = 4/5 is the point in which a minimum at (3 = starts to 
develop and defines the antispinodal line. 

• For x 7^ (y ^ 0) there exists a region, where two minima, one spherical and one deformed, coexist. This region 
is defined by the point where the (3 — minimum appears (antispinodal point) and the point where the (3 =/= 
minimum appears (spinodal point). The spinodal line is defined by the implicit equation: 

! + ^ (13, 




A, 

where A = (4-3 x + 2 (x- 1) y 2 ) 2 and B = 36 j/ 2 (ir- 1) 2 . In the SU(3) case, \ = ~V7/2, provides x ~ 0.820361. 

In the coexistence region, the critical point is defined as the situation in which both minima (spherical and 
deformed) are degenerate. At the critical point the two degenerated minima are at 0o = and (3q = a^x/^ 
(A) = J//2) and their energy is equal to zero. The critical point line can be calculated to be 

_ 4 + y 2 4 + x 2 (£,0;£0|£,0) 2 
c 5 + y 2 5 + x 2 {L,0;L0\L,0) 2 ' [ ' 

In the case of L = 2: 

~- 4+ ^ (15) 



5 + fx' 



2 ,.2 ' 



being in the 5/7(3) limit (x = -V7/2) x c = 9/11. 



• According to the previous analysis, for ^ ()/ / 0) there appears a first-order phase transition, while for 
X = 0(y = 0), there is an isolated point of second-order phase transition at x = 4/5. In this last case, 
antispinodal, spinodal and critical points collapse in a single point. 

The substitution of (3o(\o) in the Hamiltonian i|llfl implies that the term of order N 1 / 2 vanishes because it is 
proportional to the derivative of Hi with respect to A. More precisely, one has that dHi/dX = 2Hi/ 2 - The first 
quantum corrections comes from the N° term which is a simple quadratic form in the c boson operators. It cans 
thus be diagonalized through a Bogoliubov transformation. This transformation depends on the phase, spherical or 
deformed and in the next subsections both will be treated separately. 
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A. Bogoliubov transformation in the spherical phase 



In the spherical phase (3 = (A M = for all /i) and x > 4/5. In this case, the Hamiltonian reads as, 



H = (3x - 2)n c + (x- I) 



'2L+1) + (P} +P r 



0(l/N), 



(16) 



which is straightforwardly diagonalized via a Bogoliubov transformation 

4 = + vJu 



(17) 



where the coefficients verify u 



= 1, with Uu, = 



and w M = U-jn. The phases of the coefficients are chosen so 



as to minimize the mean field energy, leading to 



H 



2L + 1 



x + E(x) 1 ^ 2 ] + n£ E{x) 1/2 + 0(1/N), 



(18) 



where we have introduced E(x) = x(5x — 4) and is the number operator for £ bosons. Note that in the spherical 
phase the mean-field energy is equal to zero. In this phase, which is only defined for 4/5 < x < 1, the spectrum is, 
at this order, independent of \ (y) an d fias a trivial dependence on L. As shown in Ref. |33| for L = 0, one has to 
diagonalize H at next order (1/N) to see the role played by this parameter. 

In this phase there exists a (2L+1) times degenerated phonon (5 in the IBM case), £. The Hamiltonian is completely 
harmonic and therefore the two-phonon excitation energy is exactly twice the one-phonon excitation energy. 

Another observable of interest that can be calculated easily is the number of L bosons in each state. For the 
calculation of such observable the Hellmann-Feynman theorem can be used. It establishes that the derivative of the 
eigenvalue of a given operator, e.g. the Hamiltonian, is equal to the expectation value of the derivative of this operator 
with the corresponding eigenfunction. This leads to: 



d_ 

00 



[1 + 6){H) 



(19) 



where 9 = y^j- In this case, the contribution from the mean field is zero and the first non vanishing contribution 
comes from the term proportional to 7V° in the energy. Therefore, 



2L + 1 



3x - 2 



1 



E(x) 
"3x-2 

P 



E(x) 



hO(l/N), 
■OQ./N), 



(20) 
(21) 



where (rii) gs stands for the expectation value of the number of L bosons in the ground state, p is the number of 
excited £ bosons and {n^jp^ stands for the expectation value of the number of L bosons in the state with p excited £ 
bosons. The N° correction is singular at x = 4/5 as already noted in similar models |34ll35|. 

Note that here, we have chosen [3q as an order parameter but, one may have taken {ni) g s equivalently. Indeed, in 
the thermodynamical limit, this quantity is only nonvanishing in the deformed phase as we shall now see. 



B. Bogoliubov transformation in the deformed phase 

In the deformed phase where /?o 7^ (Ao ^ 0), the situation is more complicated and strongly depends on I. In 
the following, we will discuss the two cases L = 0, 2 separately but we underline that the form (|llfl of the expanded 
Hamiltonian allows for the study of arbitrary L. 



1. The case L = 



The case L = has recently attracted much attention because, at the mean field level, it reproduces exactly the 
IBM phase diagram (although, of course, it does not include K = 2 excitations) . In Ref. [33| , we have computed the 
finite-size corrections up to 1/7V order in the spherical phase. Here, we shall now treat the deformed phase at order 



(1/N)°. At this order, for L = 0, the Hamiltonian is easily diagonalized via a Bogoliubov transformation over the c 
scalar boson and one gets, 



H = E(x,y,p )- 



1 



-x+(7x-8)ffi+2{x-l)yp {-2-yf3o+2P 2 o) 



$(s,y,/3o) 1/2 



+n^(x,y,f3 n ) 1 / 2 +0(l/N) 7 
(22) 



2(l + /9 2 ) 
where E(x,y,f3o) is given by Eq. 

x - (3s - 4)/3 2 + 2(x - l) V p {2 + y/3 - /3 2 )| [5a - 4 - (19s - 20)/3 2 + 2{x - l)y/? (6 + 3y/3 - 7/3 2 - 
$(x,y,/3 ) = 1 



2^2 



(1 + $. 

(23) 

and y — X- I n this case, one has a single phonon excitation with K = 0. For /3q = 0, one recovers expression Ijl8|) 
setting L = 0. 

Regarding the expectation value for the number of L bosons, it can be calculated as before through the Hellmann- 
Feynman theorem (see Eq. dl). Note that in the deformed phase there is a contribution proportional to N coming 
from the mean field energy. More precisely, one has 



4 



d 



i 



-x + (7x- 8)/3 2 



+ 2(x-l)y(3 (-2-yPo + 2f3 2 ) 



^(x,y,(3 y/ 2 
2(1 -x) 



(n L )pt = (n L ) gs +p(l-x) — I ^— j 



(24) 
(25) 



where we have used the same notation as in Eqs. (|20I21|) and p is the number of excited £ bosons. 



2. The case L = 2 



In this section we will focus on the IBM case, i.e. L = 2. For arbitrary L ^ 0, the Hamiltonian must be 
diagonalized for each value of \x separately. Indeed, one has 



H o — C + H M , 

/i=— 2 



(26) 



where C is a constant and = H—^. As can be seen in Eq. depends not only on \x but also on the angular 

momentum L via the Clebsch-Gordan coefficients ctj^J,. 

We diagonalize separately the modes fi = 0, fx = ±1, and [i = ±2 which correspond to the j3 phonon (K = 0), a 
Goldstone phonon (K — 1 two- fold degenerate), and the 7 phonon (K = 2 two- fold degenerate), respectively. After 
the diagonalization via a Bogoliubov transformation, the full diagonal Hamiltonian in the deformed phase reads, 



H = E(x,y,p ) + 



1 



2(1 + /? 2 ) 



with 

®±2(x,y,p ) 



-5x+(19 a ;-24)/3 2 + 12(x-l) 2/ /3 3 ]+ £ («,A) + Op), 

/x=-2 

(27) 



[a - (3a: - 4)/3 2 + 2(x - l)j/A(2 + yA - Po)] [to - 4 - (19a: - 20)A 2 + 2(x - l)j/A(6 + 3yA - 7/3 2 - At)] 



[as - (3a; - 4) A 2 + (z - l)y#>(2 + j/A - 2/3 2 )] [to - 4 - (3z - 4)A> + 2{x - l)y/3 {3 + y/3 - A?) 



;i+a 2 ) 



[a; - (3a; - 4) A 2 - 2(z - l)y#>(2 + yA + Po)] [to - 4 - (3x - 4) A 2 + 2(s - l)yA( - 6 + yA> - A 2 ) 



(28) 



(29) 



(30) 
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where y — —^/2/7x- For /3q = 0, the symmetry between modes is restored [$^(x, y, 0) = c(x)] and one recovers the 
expression 

(HE) with L — 2. 

For /3q 7^ 0, the phonon excitations depend on /i. The excitation for fi = Q bosons, which corresponds to /3 bosons, 
is the same as in the L = case, namely $o(a ; ) 2A A)) — ^( x iVi M- I n addition, the excitation energy for ji = ±1 
modes vanishes since for /3q ^ 0, one has <&±i(x,y, flo) oc aE ^y^ | = 0. This is in agreement with the fact that 
the /i = ±1 excitation corresponds to a rotation of the ground state, i.e. to a Goldstone phonon. Finally, the \i = ±2 
excitation corresponds to a 7 excitation, which is two-fold degenerate. 

For the calculation of the expectation value for the number of d bosons we proceed as in the L = case and get: 

<•">» = "ttm + " " * )2 i ( 2d + m - x) t - fa + " te - 24)/) ° 

+ i2 (l -ih»as] +E 2i |^ 

/i=-2 \ ' 

(n L ) P ^ = (n L ) gs +p(l - x) — I — ] , (32) 

where, again, the same notation as in Eqs. I|20I21|) is used and p corresponds to the number of £ M excited bosons (/? 
or 7 bosons in the L = 2 case). 

IV. NUMERICAL RESULTS 

In this section we compare the analytical results obtained in previous sections with numerical calculations. Note 
that for clarity only the firsts + and 2 + states are plotted as members of the different bands. 

A. The case L = 

In this case, we perform the numerical calculations using the technique presented in reference |33l l>i(il | . It allows to 
easily deal with a large number of bosons, up to a few thousands. One can reach such a number of bosons due to the 
underline 0(5) symmetry which allows to use a seniority scheme reducing considerably the dimension of the matrices 
to be diagonalized. 

In Fig.|21we compare the analytical with the numerical excitation energies for a large number of bosons, N = 5000 
(all the following calculations for L = are performed for N = 5000) and \ = ~ V7/2. Note that single and two 
phonon excitations are equally well described. The left part of the figure corresponds to the deformed phase while 
the right part to the spherical one. In this case there appears a first order phase transition and at the critical point 
x c = 9/11 the ground state and the first excited state are degenerated, one corresponding to the spherical and the 
other to the deformed ground state. 

In Fig. 13 we repeat the same comparison for the case x = 0. In this case a second order phase transition appears. 
The energy of the first excited state becomes zero in the deformed and in the spherical phase at the critical point, 
x c = 4/5. At this point, regarding the analytical calculations, the deformed /3-excitation transforms into the spherical 
one-phonon excitation. However, concerning the numerical results, the 0^" state, identified with the (3 band in the 
deformed sector, transforms into the two-phonon excitation in the spherical sector. 

Note that although in the spherical phase the N° correction is independent on \ (y) there is a noticeable difference 
between Fig.|5]and Fig. [3] because for each x- value only the phase that corresponds to the lowest mean field energy is 
plotted. The spherical phase only becomes the most stable from x > 9/11 on for \ = — V7/2, while in the case \ = 
it is from x > 4/5 on. It should be noted too that in the deformed phase for x = there appear degenerate doublets 
of levels due to the extra parity symmetry in the Hamiltonian in this case. Thus, the (3 band is connected to two and 
three phonon excitation in the spherical phase while the 1 band is related to the four and five (not shown in Fig. 0) 
phonon excitation in the spherical phase. 

For the number of bosons, we compare the analytical formulae with the numerical results for the case of x — — v7 /2 
and x = in Fig.0] In particular we are interested in the study of the iV corrections for the ground state, therefore 
we subtract the mean field contribution, (no)™ * /N, to both, analytical and numerical results. As expected we observe 
how the N° correction improve the description of (no) gs specially near the critical point. 
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0.78 0.8 0.82 0.84 

X 

FIG. 2: (Color online) Behavior of one and two phonon excitation energies, in arbitrary units, for L = as a function of x near 
the critical point for \ — —s/7/2. Lines are the analytical results and dots are the numerical calculations. 




FIG. 3: (Color online) Same as Fig. [2] but for \ — 0- The excited phonon in the deformed region (a; < 4/5) is the equivalent 
to the f3 excitation in IBM. The lowest + corresponds to 0^. 



B. The case L = 2 

For L = 2, the numerical calculations have been carried out with an IBM code [32j which has been modified for 
allowing calculations up to N = 100 bosons. All numerical calculations for IBM presented below are performed for 
N = 100. 

For L = 2 the case \ = reduces to the L — situation already discussed. In particular, the analytical ground 
state energy is the same in both cases, although the iV° correction differs; there exists only one kind of excitation: 
the 0, while the fi = ±1 and 7 excitations become spurious Goldstone bosons. The excitation energy is equivalent 
to (J2HJ- On the exact diagonalization side, the Hamiltonian can be rewritten in terms of the generators of an 
SU(1, 1) algebra [33|,|3(| in the same way that in the L = case. Consequently, in this section we will only consider 
L = 2 with x 7^ 0- Any \ value can be analyzed but, as an illustration, here we will present results for the case 
X = — v7/2 that gives the U(5) — SU(3) leg in the Casten triangle. 

First, we plot the analytical results corresponding to one and two phonon excitations (Fig. [SJ). In the deformed 
phase, the bosons are and 7 excitations, while in the spherical phase they are spherical harmonic phonons. At the 
critical point the and the 2 bands transform into one and two phonon bands, respectively. However, the 7, 0~f, and 
7 2 bands apparently disappear when entering in the spherical phase. Indeed, it happens because and 7 excitations 
become degenerate for = 0. The spherical phonon excitation is a five degenerate excitation where the deformed 
and 7 excitations collapse together with the Goldstone boson with projection ±1 (which is at zero energy in the 
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FIG. 4: Variation of «n L=0 ) 9S - (n L=0 )™/)/iV as a function of a; near the critical point for L = 0. Full line is for the analytical 
and circles for are numerical results. Upper figure corresponds to % — an d lower to \ — ~ /2. 




FIG. 5: (Color online) Excitation energies (analytical), in arbitrary units, for one and two phonon states as a function of x for 
L = 2 and x = ~V7/2. 



deformed phase). 

In order to compare analytical and numerical results we will split the analysis in three different regions: deformed 
phase (Fig. |SJ, critical region (Fig. EJ and spherical phase (Fig.|HJ|. The harmonic character of the results is observed 
in all these plots. 



1. Deformed phase 

In the deformed phase (Fig. |SJ) one and two phonon excitations are clearly separated in energy. Note that the 
excitation energy for the 7 band is higher than the corresponding one for the f3 band, although for x = (SU(3) 
limit) they are degenerated. Also note that the 7 2 excitation carries the angular momentum projections K — 0, ±4 
which in this approach are degenerated. 

The correspondence between numerical and analytical states is as follows: (3 band is identified with 0^, 7 with 2g", 
/3 2 with O3 , /?7 with 2$ and 7 2 with O4 . Note that the state 2% belongs to the (3 band, while 2\ to the f3 2 band. 

The overall agreement between analytical and numerical results is satisfactory and improves the description given 
in |3l| for single and double phonon excitations, although in the present approach, no mixing between the different 
kind of excitations appear. 

The average number of d bosons in the deformed phase, normalized to the total number of bosons, is depicted in 
Fig. OH for one phonon states (the results for two phonon states are not presented for clarity). It can be observed a 
smooth decrease of (rid) when x is increasing as it is expected when approaching the spherical phase. 
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FIG. 6: (Color online) Excitation energies (analytical and numerical), in arbitrary units, of one and two phonon states as a 
function of x for L = 2 and x = —y/7/2 in the deformed phase. Lines correspond to analytical an dots to numerical results. 
The lowest + state corresponds to 0^ and the 2 + states are respectively 2^ and 2g . 




FIG. 7: (Color online) Excitation energies (analytical and numerical), in arbitrary units, of one and two phonon states as a 
function of x for L — 2 and \ = ~V7 /2 in the region around the critical point. The dots correspond to numerical results. The 
lowest + state corresponds to 0^ and the lowest 2 + state corresponds to 2+. 



2. Critical area 



The comparison around the critical area (Fig. |7J| becomes complicated because one and two phonon states have 
comparable energies and there appears interchange of character between states. For example at the critical point, the 
1 is at lower energy than the 7 excitation. 

Starting at x — 0.75 the correspondence between analytical and numerical states is similar to the one given in 
preceding section, but already at x = 0.8 different states interchange its character. The correspondence between 
states is presented in table Q] From this table it is clear that there exists an interchange of character between the 
states corresponding to the 7, /3 2 , Pj, and j 2 bands. 
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FIG. 8: (Color online) Excitation energies of one and two phonon states, in arbitrary units, as a function of x for L — 2 and 
X = — vf/2. Lines correspond to analytical an dots to numerical results. 
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FIG. 9: (Color online) {rid)/N as a function of x for L = 2 and x = —V7/2, in the deformed phase for one phonon states. 
Lines correspond to analytical and dots to numerical results. 

An interesting question that arises is if the interchange of character is due either to level crossing or to level repulsion. 
We have to take into account that the transition between SU(3) and U(5) is not an integrable path i.e. a complete 
set of mutually commuting Hermitian operators does not exist. This implies that crossings are forbidden and only 
repulsion is allowed. In particular, in the thermodynamical limit the repulsion becomes anti-crossing, i.e. infinity 
repulsion. In Fig. 1101 we show a zoom of one apparent crossing in Fig. between 2 + states in the region around 
x = 0.796, it is clearly seen that the levels indeed repel each other as expected. In order to illustrate this result and 
show how the two involved levels interchange their character, we present in Fig. ^2 the expectation value of the d 





x = 0.75 


x = 9/11 


x = 0.85 


p 


2 + ,2+ 


0^,2+ 




7 


2+ 


2+ 

z 4 




P 2 


n+ 2+ 

u 3 > z 4 


0^,23+ 




Pi 




2 6 + 




7 2 


o 5 + 


0+ 




1 phonon 






2f 


2 phonons 






0+ 2+ 



TABLE I: Correspondence between analytical and numerical states for three values of x: x 
critical point and 0.75 spherical phase. Only + and 2 + states are indicated explicitly. 



0.75 deformed phase, x — 9/11 
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FIG. 10: (Color online) Excitation energies (numerical) of 2g and 2\ states at the region of closest approach (see text) for 
L = 2 and x = —s/7/2. 
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FIG. 11: (Color online) (n^/JV (numerical) of 2g~ and 2\ states at the region of closest approach (see text) for L = 2 and 
X = -a/7/2. 



boson number in both states. It is clearly observed that the states interchange their properties at the point of closest 
approach. 

The average number of d bosons, normalized to the total number of bosons, in the region around the critical point 
is depicted in Fig. ^] for the ground state (left panel) and for the (3 band (right panel). One important feature is the 
discontinuity appearing at x c — 9/11 due to the existence of a first order phase transition. In the evolution of the 
f3 band, it appears a kink in the numerical results at the critical point. This behavior at the critical point has been 
already observed for other observables such as isomer shifts derivative of the ratios of 4^/2+ excitation energies 
HH or B(E2 : 4f -> 2^)/B(E2 : 2^ -> 0^) 15]. Also note that the 0^ state transforms into a two phonon band 
when passing to the spherical phase. 



3. Spherical phase 

The last region of interest is the spherical phase (Fig. |HJ. Here, there exist a five degenerated phonon excitation. 
The correspondence between the analytical and the numerical results is clear: one phonon excitation corresponds to 
the state 2+ while two phonon excitation to the state 0^" (also to 2% and states). 

The average number of d bosons, normalized to the total number of bosons, in the spherical region is depicted 
in Fig. 1131 for the one and two phonon states. The main discrepancies between numerical and analytical results, as 
expected, appear close to the critical point. Note that the structure of the states is very simple and already for x = 0.9 
the number of d bosons is fixed to 1 and 2 for one and two phonon states, respectively. 
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FIG. 12: (Color online) (rid) jN as a function of x, for L — 2 and \ ~ — v7 /2, at the critical area for one phonon states. Lines 
correspond to analytical and dots to numerical results. 
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FIG. 13: (Color online) {nd)/N as a function of x for L — 2 and x 
states. Lines correspond to analytical and dots to numerical results. 



-yf/2, in the spherical phase for one and two phonon 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied two-level boson models characterized by a lowest scalar s-boson and an excited L 
boson through a Holstein-Primakoff transformation that allows to treat explicitly order by order a TV expansion. This 
treatment shows that only the leading N term of the ground state energy is correct in a mean field (or Hartree-Bosc) 
approach. We stress that the equilibrium nuclear shape corresponding to an IBM Hamiltonian should be obtained 
only considering the leading TV term of the ground state energy. 

Depending on the value of L, models of interest in different fields can be obtained. Thus, L = is related to the 
Lipkin model first introduced in Nuclear Physics and then used in many fields, L = 1 is the vibron model of interest 
in Molecular Physics, L — 2 is the Interacting Boson Model of nuclear structure, etc. We have presented a method 
to go accurately beyond the standard mean field treatment so as to be able to compute finite size corrections to 
several spectroscopic observables. The model Hamiltonian used is a generalization for arbitrary L of the Consistent Q 
Hamiltonian in the IBM. This Hamiltonian depends on two control parameters and changes in them allow to explore 
the full model space and the corresponding phase diagram. Although the formalism is general for any L value, we 
have concentrated in the cases L = 2 (IBM) and L = 0. Both spherical and deformed phases have been studied with 
special emphasis in the gap for single and double excitations and the expectation values of the number of L bosons 
in different states. Analytic results have been validated by comparison with full numerical calculations. 
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